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1. INTRODUCTION 


1.1 Overview of Research Project 

The design of helicopter rotor blades involves not only considerations of 
strength, survivability, fatigue, and cost, but also requires that blade 
natural frequencies be significantly separated froia the fundamental aerodynamic 
forcing frequencies (e.g. Ref. 1). A proper placement of blade frequencies 
is a difficult task for several reasons. First, there are many forcing 
frequencies (at all integer-multiples of the rotor RPM) which occur at rather 
closely-spaced intervals. For example, 5/rev and 6/rev are less than 20Z 
apart. Second, the rotor RPM may vary over a significant range throughout the 
flight envelope, thus reducing even further the area of acceptable natural 
frequencies. Third, the natural modes of the rotor blade are often coupled 
because of pitch angle, blade twist, offset between the mass center and 
elastic axis, and large aerodynamic damping. These couplings complicate the 
calculation of natural frequencies. In fact, the dependence on pitch angle 
makes frequencies a function of loading condition, since loading affects 
collective pitch. Fourth, the centrifugal stiffness often dominates the 
lower modes, making it difficult to alter frequencies by simple changes in 
stiffness. 

In the early stages of the development of the helicopter, it was believed 
that helicopter vibrations could be reduced (and even eliminated) by the correct 
choice of structural coupling and mass stiffness distributions. However, it 
is easy to imagine how difficult it is to find just the proper parameters such 
that the desired natural frequencies can be obtained. The difficulties in 
placement of natural frequencies have led, in many cases, to preliminarv 
designs which ignore frequency placement. Then, after the structure is 
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"finalized" (either on paper or in a prototype blade), the frequencies are 
calculated (or measured) and final adjustments made. Reference (2) describes 
the development of the XII-17 helicopter in which a 300-lb weight was added to 
each blade in order to change the spanwise and chordwise mas3 distribution and 
thereby move the first flapwise frequency away from 3/rev. However, these 
types of alterations are detrimental to blade weight, aircraft development 
time, and blade cost. In addition, corrections usually are not satis- 
factory; and the helicopter is often left with a noticeable vibration problem. 

The state-of-the-art in helicopter technology is now to the point, 
however, that it should be possible to correctly place rotor frequencies 
during preliminary design stages. There are several reasons for this. First, 
helicopter rotor blades for both main rotors and tail rotors are now being 
fabricated from composite materials (Refs. 3 and 4). This implies that the 
designer can choose, with certain restrictions, the exact El distribution 
desired. Furthermore, the lightness of composite blades for the main rotor 
usually necessitates the addition of weight to give sufficient autorotational 
blade inertia. Thus, there is a considerable amount of freedom as to how 
this weight may be distributed. Second, the methods of structural optimization 
and parameter identification are now refined to the point where they can be 
efficiently applied to the blade structure. Some elementary techniques have 
already been used for the design of rotor fuselages (Ref. 5). It follows that 
the time is right for the use of structural optimization in helicopter blade 
design. Some work on this is already under development, and, although not 
published, some companies are already experimenting with the optimum way to 
add weight to an existing blade in order to improve vibrations. 


\ 
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The purpose of Che work discussed here Is to investigate the possibilities 
(as well as Che limitations) of tailoring blade mass and stiffness distributions 
to give an optimum blade design in terms of weight, inertia, and dynamic 
characteristics. 

The major objectives of the work are: 

1) To determine to what extent changes in mass or stiffness distribution 
can be used to place rotor frequencies at desired locations. 

2) To establish theoretical limits to the amount of frequency shift. 

3) To formulate realistic constraints on blade properties based on weight, 
mass moment of inertia, size, strength, and stability. 

4) To determine to what extent the hub loads can be minimized by 
proper choice of El distribution. 

5) To determine if the design for minimum hub loads can be approximated 
by a design for a given set of natural frequencies. 

6) To determine to what extent aerodynamic couplings might affect the 
optimum blade design. 

7) To determine the relative effectiveness of mass and stiffness distribution 
on the optimization procedure. 

8) To determine to what extent an existing blade could be optimized with 
minimal changes in blade structure. 

9) To develop several "optimum profiles" for rotor blades operating under 
various standard conditions. 

The work is to focus on configurations that are simple enough to yield 
clear, fundamental insights into the structural mechanisms but which are 
sufficiently complex to result in a realistic result for an optimum rotor blade. 
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1.2 Overview of Current Report 

This second seal-annual report serves two purposes: 1) it informs our 

sponsors and other Interested parties what we have accomplished during the 
last six months, and 2) it serves as an archive of data to which we expect 
to refer frequently in coming months, as our research proceeds. As a result 
of the latter, archival purpose, the report contains much information which 
is not necessarily new, but which needs to be recorded in an orderly fashion 
so that it may be easily accessed in the future. 

The first section of the report details our experience with the CONMIN 
optimization program applied to the problem of finding the optimal design 
of a vibrating cantilever beam (see Fig. 1). This section gives the results 
of parameter changes, the results obtained with various constraint forumulations, 
and the effect of allowing lumped weights at discrete points. Other aspects 
that are studied include the autorotational constraint and the effect of 
rotation. The principal conclusion of these investigations is summarized 
as follows: CONMIN works reasonably well for all problems we have considered 

so far. 

The second part of the report discusses some numerical aspects of the 
problem. One important aspect is the effect of the r imber of finite elements 
on both the frequency calculation and the sequence of optimal designs. Also 
of importance is the effect of errors in the eigenvalue analysis as well as 
the sensitivity of the frequency to small changes in blade dimensions 
(manufacturing tolerances) . The conclusion to be drawn from these numerical 
studies is that if reasonable care is taken, numerical difficulties are not 
algnifi<ant for the problems we have considered thus far. 
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2. EXPERIMENTS WITH THE CONMIN PROGRAM AND PARAMETERS j 

I 

. \ 

The CONMIN (Ref. 6) program was employed to minimize the weight of the j 

cantilever beam, described in the introduction. Before extensive optimization 
problems are solved, it is necessary to experiment with the program parameters \ 

and determine the best values for the particular class of problems at hand. 

The problem chosen for the numerical experiments involved minimization of 
weight subject to two frequency constraints and to side constraints on the 
thickness. The non-rotating cantilever beam undergoing flapping vibrations 
wa3 examined. Specifically, the following questions were to be answered: 

1) Do the gradients produced by analytical techniques match those 
obtained by finite differences? 

2) Is there any difference between results obtained with constraints on 

j 

the frequency in Hz. and those obtained with constraints on the square j 

of the circular natural frequency in (rad/sec) ? Can scaling of the } 

i 

constraint function improve convergence? 

3) Under what conditions do the starting values for thickness influence 
the convergence properties of the problem? Can an original, infeasible 
design (i.e., one which violates one or both of the frequency constraints) 
be expected to converge to a feasible, optimal design? 

The determination of the so-called "optimal design" in CONMIN is influenced 
by several important parameters. Changes in the parameters can change the 
final answer and can influence the 3peed of convergence. The most important 
parameters are ITRM, DELFUN, DAB FUN. Convergence is defined whenever ITRM 
consecutive iterations are encountered such that the values of DELFUN or DAB FUN 
(or both) are less than the stated values input to the program. The parameters 
and default values are defined as follows: 
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ITRM: Default Value - 3. Number of consecutive iterations to indicate 

convergence by relative or absolute changes (DELFUN or DABFUN 
respectively) . 

DELFUN: Default Value ■ 0.0001. Minimum relative change in the objective 

function to indicate convergence. DELFUN » ABS (1.0 - OBJ(J-l)/ 

0BJ(J)), where the objective functions for the current, J C ^, 

s t 

iterate and the previous, J-l , iterate are tested. 

DABFUN: Default Value » 0.001 times the initial objective function 

value. Minimum absolute change in the objective function to 
indicate convergence. DABFUN ■ ABS (OBJ(J) - OBJ(J-l)). 

Note that a practical criterion for convergence is employed by the program. 
Thus, slight differences in answers can be expected if problems are started 
from different initial points, or if different values of ITRM, DELFUN, and 
DABFUN are employed. 

The parameter CT is used to define whether or not a constraint is active. 

fch 

The exact satisfaction of the J constraint, G(J) r 0 is numerically unusable. 
Rather, a band, CT in magnitude, on each side of the exact zero is employed. 

The default value and formal definition of this parameter is: 

CT: Default Value - -0.1. The J th constraint, G(J), is considered to 

be active if CT G(J) ABS(CT). The value of CT is sequentially 
reduced in magnitude during the optimization process. 

The parameter THETA is called the "Push Off Factor", and is used by the 
programmed Method of Feasible Directions (Ref. 7) to go from one feasible 
design to another feasible, improved design. The default value is 1.0. 

Larger values of THETA are appropriate for highly non-linear constraint 
functions. Lower values are appropriate as the constraints approach linear 
functions. 






The parameter PHI controls how quickly an infeasible design will be 
moved in the direction of the feasible region. The default value ■ 5.0. 
Values of PHI above 5 should be employed if a feasible solution cannot be 
obtained. 


Finally, the value of ITMAX defines the maximum number of interations 
in the optimization process. If a solution cannot be obtained in ITMAX 
iterations, the program is terminated. The default value is 10. 

Herein, experiments with the parameters ITEM, DELFUN, DABFUN, CT, THETA, 
PHI and ITMAX were performed to enable definition of the numerical values 
that best fit the class of problem at hand. There are numerous other 
parameters within CONMIN, but in this study the default values for these 
other parameters were considered to be adequate. 

2.1 The Basic Problem: Two Frequency Constraints 

The problem to be considered is the non-rotating cantilever shown 
in Figure 1. Each 24 inch long element has a different thickness of 
flange. Letting t^ be the thickness of the flange nearest the fixed end, 
element thicknesses are numbered in order, such that the tip end has a 


thickness denoted by t^Q. Corresponding values of moment of inertia of area 
are denoted by 1^ though I^q. A material density and modulus of elasticity 
of 0.1 lbf/cu.in. and 10 x 10^ psi are used. 

It can be shown that the area, , weight, an( j 

_ — ~ ~ c t r (in^) c -i — _ c *. u .. 


moment of inertia, I. 


are the following functions of the thickness. 


A^ ■ 0.50 + 7.6 t^ 
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The thicknesses are constrained to lie within the range 0.05 inches 
to 1.25 inches (the value for which the cross section would be a solid 
rectangle) . 

A uniform cantilever with a thickness of 0.10 inches for each flange 
would have the following first two frequencies: 

■ 1.98 Hz 
f 2 - 12. A Hz 

To ensure that a feasible starting solution could be obtained with at least 

one chosen thickness, the frequencies are constrained to be within +0.2 Hz 

of the frequencies 1.98 and 12. A. Thus, 1. 8 f^ < 2.2 and 12.2 _< f 2 <. 12.6 Hz. 

The mathematical programming problem becomes: 

10 10 

Min.W - Min E W « 12 + 18. 2A E t 
i-1 1 i-1 1 

S.T. 1.8 Hz < f < 2.2 Hz 
12.2 Hz < f 2 < 12.6 Hz 

and, 0.05 <_ t^ £ 1.25 

i " 1,2, ...10 (A) 

The above problem has v m thickness decision variables. To preserve 
generality for later work with more complex sections, it was decided to 
use Moment-of-Inertia decision variables. The problem then becomes: 

10 

Min.2 - E Q(I ) 
i-1 1 

S.T. G x - 1.8 - f 1 < 0 

G 2 - f 1 - 2.2 < 0 

G 3 - 12.2 - f 2 < 0 

G. - f. - 12.6 < 0 
■ A Z — 

and, 0.83073 _< I <. 5.20833 

i » 1,2. ..10 (5) 



In the above formulation, the i C ^ thickness has been written as a 
non-linear function, Q, of the i 1 "* 1 moment of inertia 

t ± - Qdj) (6) 

The problem is formulated in terms of ten moment of inertia decision variables, 

1^ through I^q, four frequency constraint functions, and twenty side constraints. 
The objective function and frequency constraints are non-linear functions of 
the decision variables. 

2.1a Analytical Gradients vs Finite Differences 

The CONMIN program has the option to compute gradients of the objective 
function and constraints via Finite Differences. If possible, however, It is 
more efficient for the user to provide analytically-derived gradients. In 
Appendix I, the method used to obtain gradients is given in detail. By using 
both methods on the same problem, it is possible to provide checks on the 
derivation and programming of the analytical gradient method. 

Tables 1A and IB present the results from using the two options. An 
original design with constant thickness of 0.05 inches leads to the designs 
in Table 1A, whereas a starting thickness of 0.10 inches is the basis for 
Table IB. It is noted that the finite difference method leads to essentially 
the same results as obtained with analytical gradients. 

The four designs all converge to a common design defined bv: 

1) A thickness of element 1 of aobut 0.086 inches - 0.091 inches 

2) The thickness of all other elements are governed by the 
lower bound constraint (0.05 inches) 

3) The second natural frequency is governed by the lower bound 
constraint (fj ” 12.2 Hz) 

4) The first natural frequency is not actively constrained, and 
ranges from 1.94 to 2.00 Hz. 
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The sensitivity of the first (lowest) frequency to small changes in 
thickness will be considered later, along with the accuracy of the single 
precision routine used for eigenvalue extraction. 

It is instructive to make a further check on the solutions by running 
an optimization problem with only one decision variable (1^) and one 
frequency constraint (f 2 >. 12.2 Hz). All other thicknesses are held 
constant at 0.05 inches. The solution is started with t^ * 0.10 inches, 
and in 6 iterations converges to the results that follow: 

Ii(ti) ;1. 256 inchea^(0.0902 inches) 

f L :2.01 Hz 

f 2 :12.21 Hz 

Weight: 21. 85 lbf 

These results can be considered to be the values toward which all four 
runs should approach. Within the limits of numerical accuracy, all four 
cases do converge to the values obtained in the calibration run. 

Hereafter, all analyses are based on the use of analytical gradients. 
2.1b Various Forms of the Frequency Constraint Function 

The constraint relations, described earlier, have been 

G x - 1.8 - f <_ 0 

G 2 - f - 2.2 < 0 

G 3 = 12.2 - f 2 < 0 

G. = f„ - 12.6 = < 0 
4 2 — 

where, f^, f 2 are the natural frequencies in Hz. (7) 

The authors of CONMIN recommend use of constraint functions that are 
mutually of the same order of magnitude, and that property is satisfied 
by the formulations in Eqn. 7. 
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It Is instructive to see whether the use of eigenvalue constrairt3 

alter the convergence of the problem. 

.th 


Let X^ ■ i eigenvalue 


( 8 ) 


■ (i)^ (the square of the circular natural frequency) 

2 

If X is extracted in units of (rad/sec) , the alternate form 
of the constraints becomes 
G* - ( 1 . 8 ( 2tt ) ) 2 - \ x 

G * 2 - X x - ( 2 . 2 ( 2 tt )) 2 


G* - (12.2(2ir)) 2 - X 


C 4 - X 2 - (12.6(271)) 


i>3 


where, X^, X 2 are the eigenvalues in (rad/sec)' 


(7a) 


The G^ constraints in Eqn. 7a are no longer expected to be of the same order 
of magnitude and convergence difficulties may result. 

To test the convergence properties, the two forma of the constraint 
function are used on a problem with initial thickness o f 0.05 inches. As 
seen in Table 2, the formulation with eigenvalue constraints does converge, 
but the optimal solution is not as good as that obtained by using frequency 
constraints . 

When the eigenvalue constraints are used on a problem starting with 
t * 0.10 inches, a very poor optimal solution is obtained. A good optimal 
solution could be obtained by scaling the objective functions by 10 and 100. 

A* A 4 
- Gj^/10 



G*/10 


** * 

C 3 - G 3 /100 

A* A 

G, - C./100 
4 4 


(7b) 


1 

i 


-V. 
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Alternately, constraints could lead to a good optimum if the constraint 
thickness, CT, were changed from -0.1 (the default value) to -8C0.0. 

When the eigenvalue constraints are used on a problem starting with 
t ■ 1.25 Inches, no feasible design is found. Attempts to move toward a 
feasible design by changing the parameter PHI from 5 to 50 to 150 are 
unsuccessful. 

For all further studies, frequency constraints of the form shown 
in Eqn. 7 are employed. It is to be noted that if troubles are encountered, 
a more efficient form of Equation 7 can be employed to ensure objective 
function values of the same order of magnitude. 




In Case 1, the initial design violates the lower bound on while in 
Case 3, the lower bound on both frequencies is violated. Only in Case 2 
does the initial design result in a feasible initial solution. 

The information for the three runs is presented in Table 3. In all 
cases, feasible solutions were readily obtained, and eventually, the optimal 
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solution wa3 obtained. If the moot difficult case (Case 3) had been 
required to meet more rigorouo convergence criteria, a few more iterations 
would have resulted in an improved optimal design. 

Certainly, initial designs which are feasible and close to optimal 
are ideal. But it is possible to start with designs which are not feasible, 
and which are far from optimal. Unfortunately, such conclusions are problem 
dependent. For more severe frequency constraints, or for added problem 
constraints, it may be necessary to start with feasible or close-to-feasible 
designs in order to optimize. 

2. Id Parameters ITMAX, ITRM, DELFUN, DAB FUN 

All data collected are based on the following values of the CONMIN 
convergence control parameters. 

ITIIAX: Maximum number of iterations (default value ■ 10) 

Values used: 40, 80 

ITRM: Consecutive iterations for convergence (default value » 3) 

Values used: 3, 5, 8 

DELFUN: Relative change parameter (default value ■ 0.0001) 

Values used: 0.0001, 0.00005 

DABFUN: Absolute change parameter (default value - 0.001 tines the 

initial objective function) 

Values used: 0.011 (default value for t » 1.10 inches) 

0.001 (default value for t m 0.10 inches) 

0.0005 (default value for t » 0.05 inches). 

In general, the default values for ITRM, DELFUN, and DABFUN lead to 
good convergence properties. In about half of the runs, the objective 
function does not change, or just barely changes, during the last ITRM 
iterations. The other half have convergence governed by the DABFUN parameter. 


14- 


Whon larger thicknosaen are used, higher initial objective function reaulta, 
and tha dofault value for DABFUN can be larger than desired. 

To tent the adequacy of the default values, Case 2 in Tablo 3 cay be 
rerun with tighter controls. The results, iteration by iteration, are shown 
in Table 4. The tighter controlled run essentially doubles the ntcaber of 
iterations, and halves DELFUN and DAE FUN . Yot, the final results are 
essentially identical. 

The recommended parameters for such runs are: 

ITKAX - 40 

ITRM ■ 3 (default value) , or 5 
DELFUN - 0.0001 (default value) 

DABFUN ■ default value (0.05 inches <_ t £ 0.25 inches) 

• 0.0025 (t > 0.25 inches). 

The parameters stated may not be applicable either for larger problems, or 
for more severely constrained problems. The parameter will be critically 
examined at several stages of the study. 

2. le Parameters THETA, PHI 

Higher values of THETA (The Push Off Factor) are recommended for highly 
non-linear constraint functions. The default value (0 ■ 1), 10, 100, and 
700 were tried on a specific Initial design (t ■ 0.10 Inches I - 1.355 inches ). 
The results are shown In Tablo 5. Clearly, changing the parameter from the 
default value does not Improve the rapldltv of the convergence or the quality 
of the answer. 

In section 2.1b, experiments with PHI were reported. Changing HU from 
5 (the default value) to 50 and 150 did not enable an initially infeasible 
design to bo moved into the feasible design space. 


4 
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Further designs are hereafter based on default values for THETA and 

PHI. 

2.2 Placement of Frequencies (More Severe Constraints) 

The constraints defined by Eqn. 7 can be generalized as follows: 

G r £ iL- £ ii° 

C 2 - f l - f lU 1 0 
G 3 " f 2L - f 2 ± ° 


G 4 “ f 2 " f 2U - 0 


( 10 ) 


th 


where, the lower and upper bounds on the i frequency, 
in Hz, is given by f^ L and f respectively. 

Three cases are examined, as follows: 

Case 4: Cases 1, 2, 3 were constrained by a band of + 0.2 Hz 

around 2.0 Hz and 12.4 Hz. The band is now narrowed to 
+ 0.1 Hz. 

Case 5: The first two frequencies are separated, such that a 

+ 0.2 Hz band is defined around 1.7 Hz and 13.0 Hz. 

Case 6: The first two frequencies are brought closer together, such 

that a + 0.2 Hz band is defined around 2.3 Hz and 11.8 Hz. 

The results of the optimizations are shown in Table 6. Thus, within reason, 
it is possible to re-proportion initial designs (feasible or infeasible) such 
that frequencies are placed where desired during an optimization of weight. 
2.2a Further Studies of Convergence Parameters 

The recommended ITMAX, ITRM, DELFUN, and DABFUN parameters (30, 3, 
0.0001, 0.0025) are used in Case 4. The convergence of the objective 
function to 3 significant figures seems to be incomplete in the 3rd figure. 
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For Case 5, ITRM is changed from 3 to 5 and a surprisingly substantial 
reduction in objective function is obtained. Five iterations appear to be 
appropriate to maintain objective functions of 3 significant figures. For 
Case 6, another check (not shown in Table 6) wa3 made by tightening the 
parameters to 40, 6, 0.00005, 0.0001. After 40 iterations, the convergence 
criteria had not been satisfied, but probably would be in another few cycles. 
The objective function was only changed from 0.593 to 0.590, but elements 
2 and 3, originally 0.0703 and 0.0886 inches, were appreciably changed to 
0.0832 and 0.0791 inches respectively. The objective function is quite flat 
near convergence, and small changes in objective function can be accompanied 
by appreciable changes in structural configuration. The accuracy of 
eigenvalue extraction, herein done by a standard library routine in single 
precision, obviously has an effect on the defined final configuration. 

A numerical study of eigenvalue calculations is presented in section 5. 
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3. OPTIMAL DESIGNS FOR A NON- ROTATING CANTILEVER 


With the experience gained by studying two frequency constraints, it 
is possible to intelligently formulate the more difficult problems that 
follow: 

1) Added Constraints - A third frequency 

2) Added Design Variables - Non- structural lumped mass 

3) Added Constraint - The auto-rotation constraint (i.e., minimum 
mass moment of inertia about an axis normal to the beam, and 
located at the root) . 

3.1 Optimization with Three Frequency Constraints 

In addition to the constraints shown in Eqn. 10, two more must be 
added: 


where, 
in Hz, 


C 5 * f 3l - f 3 i 0 


°6 ‘ f 3 ' f 3U ± 0 


the lower and upper bounds on the 3rd frequency, 


is given by f^ and f ^ respectively 


(ID 


Three cases are examined, as reported in Tables 7A and 7B: 

Case 7: A + 0.2 Hz constraint band is placed around the first three 

frequencies of the initial design. 

Case 8: The spread between the desired values of f^ and f ^ is 

narrowed, and the desired value of f^ is decreased by 
1.0 Hz from that defined in Case 7. 

Case 9: The spread between the desired values of f^ and f^ is increased 

and the desired value of f^ is raised by 1.0 Hz from that 


defined in Case 7. 
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The more difficult optimizations (Cases 8 and 9) are repeated with different 
starting points, and (in Case 9) with different convergence criteria. Con- 
vergence to an estimated 99% accurate value of objective function is 
accomplished for runs using 1TMAX • 80, ITRM ■ 5, DELFUN ■ 0.0001, and 
DABFUN - 0.0005. Designs initiated with various constant thickness values 
converge to similar optimal designs. For example. Cases 9B and 9C both 
result in minimum thickness for elements 2, 3, 4, 5, 9, and 10. Despite 
identical objective functions, individual thicknesses for elements 6, 7, 
and 8 in Cases 9B and 9C are far from identical. Again it is noted that 
near optimum there can be appreciable changes in structural configuration 
with very minimal effect on the objective function. 

3.2 Addition of Lumped Weights 

The introduction of lumped weights (assumed to have mass, but not to 
contribute to the moment of inertia of area, I, used in defining the stiffness 
matrix) increases the number of decision variables in the optimization 
problem. For one such weight at the center of each element, there are now 
twenty decision variables: 

I^(root element) through (tip element) 

W^ (weight on root element) through W^^ (weight on tip element) 
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The problem is formulated as follows: 

10 10 

Min Z - 182.4 y Z Q(I ) + Z W 
i-1 1 i-1 1 

S.T. - f 1L - f 1 £ 0 


f l- 


f iu±° 


2L 


- f £ £ ° 


g 5 - f 3L - f 3 < 0 

G 6 " f 3 - f 3U - 0 


and, 0.83073 £ I ± £ 5.20833 

i-l,2,...10 

0 < W ± < 100 

i**l,2, . . .10 


( 12 ) 


The newly defined symbols in Eqn. 12 are as follows: 

y “ element density, lbm/cu.in. 

til 

«• weight of lumped mass at i element 
center , lbm 

Very minor modifications are necessary to include the analytical 
gradients of the objective function with respect to the lumped weights 
and of the constraints with respect to the lumped weights: 






( 1 * 1 , 2 , . . . 10 ) 


( 13 ) 



where, Z is the objective function 

K is the stiffness matrix 
M is the mass matrix 
I is the identity (unity) matrix 

2 

g is the acceleration of gravity, 386.4 inches/sec 
The two matrix derivatives in Eqn. 13 are used to define the partial deriv- 
ative of the eigenvalue with respect to the decision variable, OX^/SW^). 

The remainder of the operations are outlined in Appendix I. Thus, the 
optimization con still be based on analytical gradients. 

Three cases have been investigated. The results are shown in 
Tables 8A and 8B. 

Case 10: A previous design with t - 0.10 (Case 7) is modified by 

reduction of the density from 0.10 Ibf/cu.in. to 0.05 Ibf/cu.in. , 
and by replacement of lost weight by 1.512 lbf lumps at each 
element center. There is more freedom to choose the decision 
variables in Case 10 (i.e., the non-st rue tural mass can be used 
efficiently to move frequencies). The final result is a 
constant-section beam of minimum thickness with lumped weights 
as shown in the Table. The optimal weight is only about two- 
thirds that of the optimal weight for Case 7. 

Case 11: Two runs were made with only one difference in initial design. 

In 11A, 1 lbf lumped weights were used at each node, and in 
11B lumped weights were not used. The final optimal designs 
were almost identical in all respects. Since all of the initial 
frequencies had to be raised, it was most efficient to remove 
all of the lumped weights of Case 11A. 








Case 12: The Initial design of case 11A was used with a compressed 

range of frequency constraints (from 1.8-35 to 2.6-32 Hz). 

The optimal design involved a combination of thicknesses 
greater than the minimum, and lumped weights. The tip 40Z of 
the beam was made of minimum thickness elements without 
lumped weights. 

Convergence for all the runs was excellent. Note that since the objective 

function is now weight (with numerical values of 20 to 30 lbs) previouslv 

used values of DABFUN are not appropriate. Herein, DABFUN was raised to 

0.001. The other parameters were kept the same as before. 

3.3 Addition of Auto-Rotational Constraint 

This constraint is intended to be applied to rotating systems. However, 

the constraint is added here as the next step in developing the larger 

problem to be considered. Denoting the minimum mass moment of inertia about 

a vertical axis through the root of the beam as I , , and the actual moment 

min 

of inertia as I , it is required that I > I . . Rearranging the information 
m m — nin 

into a more usable form, the seventh constraint, to be added to Eqn. 12 
becomes: 

G 7 - 1 - (I /I , ) < 0 (14) 

/ m min — 

The gradients of with respect to the decision variables are analytically 
obtained, as shown in Appendix II. 

Twc optimization runs are made. In the first run, Case 10 is re-run 

2 

with a I . value of 500 lbf inches sec . This value is deliberately chosen 
min 

to be low, such that the constraint remains inactive throughout the run. 

As expected, the results remained identical to that in Case 10. The optimal 

2 

design has an I value of 807.7 lbf inches sec . The next case is reported 






~:l ■aV.'j .... 
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2 

in Table 9. A value of I . • 1100 lbf Inches oec is demanded, such Chat 

tain 

the auto-rotational constraint is active during the optimization. For com- 
parison, the new Case 13 and the previous Case 10 (without auto-rotational 
constraint) are included in Table 9. Although satisfaction of the formal 
convergence cri'.eria is not met in 80 iterations, the optimization is close 
to being finished. (DF.LFUN and DABFUN <_ 0.0004 and £ 0.005, respectively, 
for the last five iterations). As expected, the tip element is thickened 
and the lumped weight increased, since that is the most efficient way to 
satisfy the auto-rotational constraint. At optimum the auto-rotational 
constraint was, for all practical purposes, one of the active constraints. 
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4. OPTIMAL DESIGNS FOR A ROTATING CANTILEVER BEAM 

If the cantilever bean is rotating in a horizontal plane, centrifugal 
forces are created which stiffen the system and increase the natural 
frequencies in vertical vibration. The first natural frequency has as a 
lower limit equal to the speed of rotation. Thus, for high speed of 
rotation, placement of frequencies must be done with due consideration of 
the centrifugal effects. These trends are shown in Table 10. For 300 RPM, 
one bound on the fundamental frequency is 5 Hz. Despite vast ranges in 
thickness and lumped weights, f^ ranges only between 5.1 and 5.7 Hz for a 
uniform cantilever. Frequency placement is much less dependent on stiffness 
and mass distribution than for a non-rotating beam. The calculations in 
Table 10 are based on Ref. 8. 

The optimization formulation remains the same as shown in Eqn. 12 
and supplemented by Eqn. 14. The major modification required involves the 
contribution of element tension to the element stiffness matrices. The 
added contribution also causes a modification of the frequency constraint 
gradients. Some details of the derivations are shown in Appendix III. All 
eigenvalue calculations are based on double precison routines. 

In Table 11A, results for three optimizations are shown. The single 
difference in input is the speed of rotation. 

Case 14: This run is a repeat of the single precision run of Case 10. 

The non-rotating beam solutions are almost Identical. The 
major difference is a slight re-errangement of the lumped 


weights. 
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Case 15: The speed of rotation Is the low value of 30 RPM. There is 

no difficulty in placing the frequencies in the same range 
as required for the non-rotating beam. The optimization 
results in a slight increase in weight over that for the 
previous case. 

Case 16: With a speed of 100 RPH, the problem does not fully converge 

in 80 iterations. The objective function appears to be 
accurate to two (rather than the requested three) decimal 
places. The tip element and tip weight are now relatively 
large. The extra mass is needed to lower the frequency and 
counteract the effect of the high speed. 

With speeds of 300 RPM, the requested frequency placement 
must be modified. Also, the convergence criteria is relaxed, 
as shown in Table 11B. 

Case 17: This 300 RPM run converges in 39 iterations. The convergence 

criteria is satisfied by the DABFUN requirement (3 consecutive 
iterations with absolute change in objective function <_ 0.01). 
Again, a large tip mass is needed, but that may be required 
to satisfy the more demanding auto-rotational constraint. 

The use of lumped mass appears to be the more efficient way 
of controlling the fequency placement. In particular, masses 
are placed at the tip and near the zero points of the 2nd 
and 3rd mode shapes. In contrast, changes in thickness modify 
both stiffness and mass, and are less effective in "arturbing 




the frequencies 
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5. NUMERICAL ASPECTS OF THE PROBLEM i 

j 

The results reported In Sections 2, 3 and 4 are encouraging, since 
they demonstrate that, given a mathematical representation of a cantilever 

i 

beam under various conditions, the CONMIN program can produce an improved j 

design. That is, at least for the problems considered thus far, numerical 

optimization is possible. Clearly, however, the design found through the •' 

use of CONMIN will be of no use if the mathematical representation of the 

structure is at fault. Thus, in the present section, several numerical 

aspects of the accuracy of the analysis model are considered. 

5.1 Convergence with Increasing Number of Elements 

One aspect of the mathematical representation of a structure with the 
use of finite elements is the question of how many elements are required 
to obtain an acceptable accuracy. For the present optimal design studies, 
this question takes two forms: 1) are enough elements U3ed to predict the 

frequencies accurately, and 2) are enough elements used to describe the 
optimal design (that is, will essentially the same optimal design result 
if the mesh is refined)? The first question is addressed in Section 5.1a; the 
second in Section 5.1b. 

5.1a Convergence of Frequency 

To study convergence of frequency with increasing number of elements, 
a vibrating cantilever beam is considered. The beam is non-rotating. The 
analytical solution for the frequency in cps is known to be (Ref. 9) 


•y 



2ttL“ 


( 15 ) 



The results of the finite element analysis for various numbers of 

A 

elements are shown in Table 12. It can be seen that as few as six elements 'j 

.1 

gives a good approximation (less than one percent error) for the first ' 

.j 

frequency. As would be expected, the approximations for the second frequency 
are not as accurate, but the error is only about two percent when ten i.lements 
are used. In general, these re3ult3 indicate that the choice of ten elements 
in the optimization studies described in Sections 2,3 and A is justified. 

5.1b Convergence of Optimal Design 

To study how the optimal design changes as the number of elements 
increases, a cantilever bean with "N" elements and with lumped weights added 
at the nodes but otherwise similar to the beam of Section 5.1a and Figure 1 
Is considered. The density, and the constraints on the natural frequencies, 
lumped weights, and moments of inertia are 


- -w 
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Y - 0.05 lbf /in 
1.0 < ^ < 1.3 (Hz) 

10.0 < f 2 < 11.2 (Hz) 

0.0 £ \ 1 100.0 (lbf) 

0.83073 < I jL £ 5.2083 (in 4 ). 


The initial design is 


I ± - 1.3555 

W ± - 15.120/N; I - 1,2 N 


which implies an initial value of OBJ (- the objective function » the 
total weight) of 30.22A9 lbf. 

Results of the study are shown in Figures 3-8. In all cases, the 
active frequency constraints were found to be 

f x - 1.3 (Hz) 

f 2 - 11.2 (Hz) 

Figure 3 demonstrates, as one would expect, that the optimum weight 
does in fact decrease as more elements are added to the mesh. The change 
in optimum weight is quite small (note that the scale of the vertical 
axis begins at 20.0). 

Figures A and 5 show the variation of the lumped weight and the moment 
of inertia (of the cross-sectional area) at the free end versus the total 
number of elements in the mesh. It appears that these quantities do not 
converge. This result can be explained, however, by referring to Figure 6, 
in which the upper curve represents the total weight at the free end. 
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(H(N) ia the non-structural, or, lumped weight; M(N) la the structural 
weight associated with the mass distributed throughout element "N"). It 
can be seen from the figure that the total weight appears to converge 
smoothly as the mesh ia refined. The explanation for the apparent non- 
convergence shown in figures 4 and 5 and the convergence shown in the top 
curve of Figure 6 is that the "structural weight" at the free end of the 
cantilever is not really structural, since there is no portion of the beam 
beyond the free end which needs to be supported. Thus, the optimization 
routine is indifferent to whether structural or non-structural weight is 
present at the free end - the only thing that counts is the total weight 
at that end. 

Figure 6 also shows the variation of the lumped weight slightly beyond 
the middle of the beam. (All optimal designs had non-zero lumped weights 
there and at the free end of the beam.) The weight can be seen to decrease 
smoothly as the mesh is refined, although no asymptote appears present. 

A possible explanation for this behavior is that as the mesh is refined, the 
weight in the middle is being placed more efficiently - and thus less ia 
needed. 

The various sketches in Figure 7 show the distribution of mass and 
stiffness along the beam for increasing numbers of elements. It is 
interesting t observe that the optimization routine finds it most efficient 
to meet the contraints on frequency by varying the lumped weight rather than 
bv varying the stiffness (moment of inertia), since this latter quantity is 
at its lower bound everywhere except near the end of the beam. 

As was pointed out previously in reference to Figures 4 and 5, the 
optimization algorithm appears to treat the structural and non-structural 
mass at the end of the beam as interchangeable. To test this hvpothesis 
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further, the optimal design problem statement was altered slightly by 
decreasing the upper bound constraint on tho moment of Inertia from 
5.2083 to 2.0. The resulting optimum design is shown in Figure 8, and 
should be compared with the design (for N ■ 10) shown in Figure 7. Note 
that the constraint on the moment of Inertia for element 10 is not active 
in the optimal design of figure 8 (the constraint was active during the 
CONMIN iterations leading to this optimal design). Thus, the effect of 
the constraint is to lead the optimization algorithm along a different 
path than that followed when the constraint value was 5.2083. The design 
found, however, has about the same total weight at the free end (■ 9.9705 lbf) 
as the previous ten-element optimum (- 9.9222 lbf). This result confirms the 
hypothesis that CONMIN increases the moment of inertia at the free end only 
as a means of increasing the mass there. Once that option is closed (that 
is, the upper bound constraint is reduced to a value of 2.0), CONMIN simply 
increases Che lumped weight at the bean tip. This finding suggests that, 
in future optimization studies, a tight constraint be imposed on the moment 
of inertia at the free end, since little structural capability is needed 
there, and necessary end mass can be adequately represented by the lumped 
weight design variable. 

5.2 Accuracy of Eigenvalue Calculations 

Among finite element analysts, the problem of the static analysis 
of a cantilever beam subjected to an end load is notorious for being 
numerically ill-conditioned. This ill-conditioning also becomes apparent 
in Che eigenvalue calculations associated with the present optimal design 
studies. Table 13 illustrates the magnitude of the errors arising in the 
eigenvalue calculations for a non-rotating cantilever like that of Figure 1, 
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3 

with thickness 0.1 inches, density 0.1 lbf/in , and no lumped maos. Column 1 
in the table contains the first twenty frequencies of the bean, which were 
obtained by a double-precision version of a code based on the Sturm-sequence 
method with inverse iteration. The eigenvalue is defined as the square of ..he 
frequency 0 ) (in rad/sec) in the equation 

([K] - u 2 [H])(Y} - 0 . (17) 

Here, K is the stiffness matrix, M the mass matrix, and Y the eigenvector. 
Columns 2 and 3 contain the eigenvalues for the same problem, but found by 
the IBM scientific subroutine program "NR00T" (based on the Jacobi method) 
in a single-precision version (column 2) and a double-precision version 
(column 3). The difference in the first entries in columns 1 and 2 is about 
three percent. 

An alternative manner of formulating the eigenvalue problem is to 
write it as 

([M] - ^ (K ) ) {Y } - 0 . (18) 

0) 

The eigenvalue is now defined to be the reciprocal of the square of the 
circular frequency. For thi3 formulation of the problem, column 4 gives 
the frequencies found by the Sturm-sequence method, and column 5 gives the 
frequencies found by the single-precision routine "NROOT". The two methods 
now give essentially the same frequencies. The improved performance of the 
single-precision "NROOT" routine is attributable to the fact that the 
accuracy of the Jacobi algorithm is dependent on the order (in terms of size) 
in which the eignevalues are found. By contrast, the Sturm-sequence method 
is independent of the "largeness" or "smallness" of the eigenvalues. 
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The significance of these errors in the eigenvalue calculations can 
be seen by inspection of Table 14, in which are given optimal designs found 
by CONMIN using both the inaccurate eigenvalue calculation and the accurate 
eigenvalue calculation. The example corresponds to Case 6 of Table 6. The 
designs are seen to differ appreciably. 

5. 3 Sensitivity of Frequency to Small Changes in Thickness 

Because rotor blades can be manufactured only to within certain 
dimensional tolerances, the question naturally arises as to the sensitivity 
of the natural frequencies of the blade to small changes in blade dimensions. 
Clearly, if small changes in blade dimensions produce large changes in 
natural frequencies, then designing a theoretically optimum blade is futile: 
the small variations in blade dimensions introduced during manufacturing would 
destroy the optimally designed vibratory behavior. To study this question, 
the data of Table 15 were generated for the cantilever beam described in 
Section 4 (Case 14 of Table 11A) . The first column in the table contains 
the derivative of the fundamental frequency with respect to the thickness 
of the first, second, third,..., and tenth (free end) element. If 0.01 inches 
is taken as a representative manufacturing tolerance, then the data in column 1 
show that the maximum corresponding change in the first fundamental frequency 
is only about five percent (=0.01 + 4.791). Another aspect of this sensi- 
tivity study is illustrated by the data in columns 2 and 3. Column 2 
contains the derivative of the fundamental frequency with respect to the 
weight of the individual elements. Since 


and 


3f L 3f : 3t i 

35^ 3t^ 35^ 

3fn. 

3^- 9 - 12 

(t^ is the thickness; see Fig. 1) , 


( 19 ) 



-32- 


the entries in column 2 are derived by dividing the entries of column 1 by 
9.12. Column 3 contains the derivative of the fundamental frequency with 
respect to the lumped weight of the individual elements. If 0.01 lbs is 
taken as a representative manufacturing tolerance, then the data in columns 

2 and 3 show that the maximum corresponding change in the first fundamental 
frequency is about half a percent (-0.01 + 0.525). 

It is also interesting to compare corresponding entries in columns 2 and 

3 and to observe that they disagree near the blade root, with the discrepancy 
diminishing appreciably as the free end is approached. The explanation for 
this behavior lies in the fact that as the distributed weight of an element 

is increased, the cross-sectional area must of course also increase. Thus the 
second column in the table represents changes in stiffness as well as in mass. 
It follows that a unit increase in distributed weight (as represented by 
column 2) will produce a larger increase in frequency than is caused by a 
unit decrease of lumped weight (column 3), 3ince, qualitatively speaking, 
stiffness appears in the numerator and mass in the denominator of the 
frequency expression, yk/M. Thus 



The inequality is large near the root, because increasing element 
stiffness has a large effect on structure stiffness. At the tip, the 
structure stiffness is barely changed (as witnessed by the last two 
entries in columns 2 and 3: - 0.227 and - 0.229). 



s 
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6. CONCLUSIONS 


For the types of beam vibration problems described in this report, 
the following conclusions may be drawn. 

1. Either analytical or finite-difference gradients can be used. The 
finite-difference approach will of course require many more function 
evaluations, but it is easily implemented. Furthermore, the problems 
considered are not especially sensitive to numerical error in gradient 
calculations, therefore, we use analytic gradients; but we reserve the 
finite-difference gradients as a viable alternative should they be needed. 

2. The frequency constraints should be formulated directly in terms of 
frequencies, rather than eigenvalues (frequency squared). 

3. CONMIN can find an optimum design, after starting with an initially 
infeasible design, at least for the problems considered. 

4. The default values for CONMIN appear adequate, with the exception 
of DABFUN ITMAX, for which a value of 40-80 works well. 

5. Optimum designs can be found even for relatively tight constraints 
on frequencies. 

6. The objective function is relatively flat near the optimum, and 
appreciably different distributions of stiffness may yield essentially the 
same value of the objective function. 

7. Two and three frequency constraints can be handled. 

8. Using both stiffness and lumped mass as decision variables presents 
no special difficulties. Indeed, the admittedly limited experience gained 
thu3 far with these types of problems indicates that the optimization seems 
to proceed more rapidly (fewer iterations to obtain convergenc) if CONMIN 
can add lumped mass rather than just add mass in the form of structural mass. 
Thus the greater complexity of the problem, i.e., the increased number of 



decision variables, appears to be more than compensated by the greater 
freedom in choosing a design. 

9. Optimal designs for rotating beams subject to frequency and auto- 
rotational constraints present no difficulties. Again, the use of lumped 
mass (rather than structural mass) appears to be a more efficient way of 
controlling frequencies. 
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7. NEAR FUTURE PLANS 


At this writing, current research effort i3 concentrated on developing 
a generic cross-section sufficiently general that both the bending and 
torsional stiffnesses of currently existing blades can be matched. When a 
suitable generic section has been attained, it will be used in studies of 
the optimization of a non-rotating cantilever experiencing coupled flap, lag, 
and torsional vibrations. Some preliminary studies of torsional and in-plane 
vibration have already been done. 

Preliminary work has also been done on the feasibility of using an 
objective function involving the sum of the squares of frequencies, rather 
than the weight. The results look promising, and this approach will be 
pursued, especially if difficulties arise with the weight-objective function 
approach, as more complex problems are analyzed. 

Another topic of research in the immediate future will be the inclusion 
of a stress constraint in the problem of the optimization oi a non-rotating 
cantilever. Preliminary exmaination of present optimum designs show no 
particular stress problems; but, in principle, the stress constraint should 
be necessary to prevent elimination of too much material from blade designs. 
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APPENDIX I 

GRADIENTS OF OBJECTIVE FUNCTION AND CONSTRAINTS 


When analytical gradients are utilized in CONM1N, the following 
derivations are necessary. 

Let 2 ■ objective function 

gj ■ j th frequency constraint function 
t^ ■ i C ^ element thickness 


T i" 1 


th 


element moment of inertia 


. th 


1) To find 


X^ m i eigenvalue of the vibration problem 
3Z 


3Ij 


we employ the chain rule 


3Z 3Z 3t j 3Z , 3I i 

3l i “ 3t f X 3I ± " 3t '3t 

10 

Since Z - I t, 
i-1 1 


( 1 ) 

( 2 ) 


and, I< - ~ + l 4 {285t i ~ 228t i + 60. 8t 3 } (3) 

37 1 9 

Then, ~ - l{~(285 - 456t i + 182. 4t x } (4) 


Furthermore, t^ ■ ^ 

37 i 2 — 1 

Finallv , ^ (285 - 456q (I } + 182.4Q (I } } 1 


3g 


2) To find -rr* 1 - we tlso employ the chain rule 

3i k 


Hm 3g, 3f 3X t 

■ ■ - n ■ * 3 V * - v ' ■ 

3i k 3f i 3A 1 3i k 


(5) 




( 6 ) 





t ~-tt* 


- 39 - 


ORiG!i\‘AL i-ags ij 
OF POOR QUALITY 


For constraints chat are linear functions of frequency, 

?8i 

- +1, -1, or 0 


( 7 ) 


Since 


f i ’ fiT X i 


1/2 


3f l 

"^1 




( 8 ) 


oil 


• is obtained aa a function of 

the eigenvectors, mass matrix, and 
. stiffness matrix of the problem ^ 


( 1 ) 


Fox, 


R. 


Optimization Methods for Engineering Design, Add ‘son-Wen ley , 


Read ing 


1’A, 1971. 
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APPENDIX II 


ORIGINAL IVtC.: U 
OF POOR QUALITY 


GRADIENTS OF THE AUTO-ROTATATIONAL CONSTRAINT 


For a system of uniform elements and lumped weights, the mass moment 
of inertia, I , can be written as: 

tfl 


I " 

m 


E {/r 2 dm, } + 


all elements 

TTj 


, . • E (r 2 m. } 

all weights 


W, 


- E (Z 1 pA i r 2 dr} + E { ^ (rj + r 2 ) 
r i 


( 1 ) 


1 ° , 3 3 U W i , 2 2. 

" ± l ± { " r L )} + 2i (r U + r L> 


where. 


p - mass density 


A^ - area of ith element - 0.5 + 7.6 t^ 


- weight of ith lump 


g - gravitational constant 

r tI ,r T - radius from beam root to the upoer and lower end 
of the element respectively. (Note that half of 
the lumped weight has been placed at the lower end 
and upper end of the element). 


3g 7 

1) To find -rr— wc emplov the chain rule 
3I i 


31. 


i 91 

I , '31 . ' 

min i 


( 2 ) 


but. 


31 31 3A, 3t, 

m m I 1 

3! i " 3A i x 3c 1 x ai ± 


{ ^ (r 3 - r 3 ) } x (7.6) x ( 


285-456t +182. 4t 2 _1 

i L_i 


24 


(3) 
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3gy 

2) To find -rjj— , direct differentiation can be employed 


!h 

aw. 


2 2 
r U + r L 
2g 
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( 4 ) 
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APPENDIX III 

MODIFICATIONS TO ACCOMMODATE ELEMENTS IK TENSION 

For an element under constant tension, T, the potential energy is 
given by 

U " 2 ; o ( it )2 dX (1) 

where W represents the vertical deflection coordinate in the Y direction, 
as shown in Fig. 2. 

From Ref. 2, the displacement function for W is given by: 

V, , 2 a ^2 S 11 

W(X) - -- (2X - 3LX + L J ) + -f(X J - 2LX - L X) 

L L 


+ — (3LX 2 - 2X 3 ) + -|(X 3 - LX 2 ) (2) 

L J L 


where, the nodal translations are and V^, and the nodal rotations are 


V_ and V^. 

Substituting the derivatives of (2) into (1), performing the integration 

1. T 

over the length, and comparing the results to U ■ -^{v} (K^] leads to 

the identification of the added contribution to the element stiffness 
matrix, [K^.J 



1.2/L 

0.1 

-1.2/L 

0.1 

lK m l - T 


2L/15 

-0.1 

-L/30 

Symm. 


1.2/L 

-0.1 


2L/15 


(3) 


The element tensions are easily defined by computing the centrigual force 
for each lumped mass, and accumulating the r.or.al tension from tip to root of 


( 2) 

Peters, I). A, Ko, T. , Korn, A., and Rossow, M.P., First Semi-Annual Status 
Report on Design of Helicopter Rotor Blades for Optimum Dynamic Character- 
istics ,NASA-Langley ('rant No. NAO-1-250, Sept. 15, 1 Q 82. 
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the beam. To keep a constant tension in each element, the distributed and 
lumped weights were divided by two and placed at the element nodes. 

To compute analytical gradients, the partial derivatives of 
with respect to the decision variables are needed. 

*( T ) 

1) To find v 1 v ' for any element, the chain rule is used: 

K 3K 3K 3T at. 

31 " 3T X 3t X 31 

The first term represents the matrix elements of Eqn. 3. 

The middle term can be obtained by writing the element tension as the 
2 

cumulated sum of m^Q r^ terms, where, ft - rotational speed in rad/sec. Since 
each mass is a function of thickness, the derivative of tension with respect 
to thickness can be obtained. 

The final term has been defined in Eqn. 3 of Appendix II. 


2) To find 


for any element, we use: 


9k (t) 9K (T) 3T 

3W r 3T X 3W R 


The second term can be obtained by writing the element tensions caused 
by the lumped masses, and then taking the appropriate partial derivatives. 



JO. Tables 


ORIGINAL FA 3 C >* 
OF POOR QUALITY 

Initial Values 
Constant t inches: C.05 

4 

Constant I inches : 0.831 

(a) Objective Function inches: 0.500 

(b) Weight lbf: 21.12 

f^Hz: 1.86 

£ 2 | Hz: 11.6 


(c) 

Analytical Gradients 

15 

0.540 

21.85 

2.00 

12.20 (d) 

1.250(o.0896) 


Optimal Solution 

No. of Iterations 
Objective Function Inches 
Weight lbf 
f : Hz 
f 2 Hz 

▼ 1 l 4, inches. 
I^inches (t^ ) 


(c) 

Finite Difference 
Gradients 
21 

0.539 

21.83 

1.99 

12.20 

1.253(0.0899) 


(Elements 2 through 10 were essentially 
all at the lower bound value of side 
constraint on thickness: 

4 

t - 0.05 inches, I - 0.831 inches ) 


Notes: 

(a) Objective Function } inches 


10 
Z t 
i-1 


i 


10 

(b) Weight lbf - 12 + 18.24 Z t. 

1-1 1 


(c) Convergence criteria were the same, except that ITR.V was raised 
from 3 to 3 cycles when finite differences were used. 

(d) Both solutions have an active frequency constraint: f is at 

its lower bound value. 


Table L\: Comoarison of Analytical Gradient vs Finite Difference Solution 
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0RIG1NAL FAC'!- TJ 
OF POOR QUALITY 

Initial Values 


Constant t inches : 0.10 

A 

Constant I inches : 1.355 

(a) Objective Function inches : 1.000 

(b) Weight lbf : 30.24 

Hz ; 1.98 

f 2 Hz : 12.4 


(c) 


Optimal Solution 


Analytical Gradients 


Cc) 


Finite Difference Gradients 


15 

No. of Iterations 

18 

0.541 

Objective Function inches 

0.543 

21.87 

Weighty lbf 

21.90 

1.94 

fl,Hz 

1.94 

12.20 (d) 

f 2| Hz 

12.20 

1.262(0.0908) 

I 1 inches 4 (t 1 inches ) 

1.214(0.0861) 

0.331(0.0500) 

t j u 4- inches.. 
I^inches (t^ ) 

0.852(0.0519) 

0.831(0.0500) 

_ , , 4. inches. 

I^inches (t^ ) 

0.880(0.0545) 


(Elements 2-5 and 8-10 were all at the lower hound 
value of side constraint (thickness): t - 0.05 inches, 

A 

I - 0.831 inches ) 

Notes: 10 

(a) Objective function inches « I t. 

i-1 

10 

(b) Weight lbf - 12 + 18.24 E t 

i-1 1 


(c) Convergence criteria for the two runs were identical. 

(d) Both Solutions have an active frequency constraint: f 0 is at 

it lower bound value. 


Table IB: Corcnarison of Analytical Gradient vs Finite Difference Solution 



Initial Values 


Constant 

t inches : 

0.05 


Constant 

4 

I inches : 

0.831 


(a) Objective function inches : 

0.500 


(b) Weight lbf : 

21.12 


• f j Kz 

• 

• 

1.86 


f 2 Hz 

• 

11.6 


Optiml Solution 

Eigenvalue Constraints 

(c) 

Frequency Constraints 

12 

No. of Iterations 


21 

0.551 

Objective function Inches 

0.539 

22.05 

Weight lbf 


21.83 

1.94 

f 1( Hz 


1.99 

12.20 (d) 

f 2 ,Hz 


12.20 

1.150(0.0799) 

I^inches^t^ 1163 ) 


1.253(0.0399) 

0.846(0.0514) 

t . i- 4. inches. 
1^, inches (t^ ) 


0.831(0.0500) 

0.910(0.0573) 

▼ j inches. 

I, inches (t, ) 

O 1 0 


0.831(0.0500) 

0.922(0.0584) 

_ . , 4 , inches. 

I^ ( inches (t^ ) 


0.831(0.0500) 

0.875(0.0541) 

_ . . 4. inches. 

I 8/ inches (tg ) 


0.831(0.0500) 


(Elenent3 2-4, 9, 10 were all at the lower bound value of 

4 

aide constraint on thickness: t - 0.05 inches, 1 ■ 0.831 inches ) 


Notes: 

(a) Objective function } inches 


10 
Z t 
i-1 


i 


10 

(b) Weight lbf - 12 + 18.24 Z t. 

i-1 1 


(c) Convergence criteria for the two runs were identical. Gradients were 
computed by analytical techniques for the eigenvalue constrained run, 
whereas finite difference techniques were used for the frequency 
constrained run. 

(vi) Both solutions have an active C requencv constraint: f 2 is at its 

lower bound value. 


Table 2: Comparison of Eigenvalue and Frequency Constraints 
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Lf) 


Initial Values 

Case 1 

Case 2 

Case 3 

Constant t, inches 

0.05 

0.10 

1.101 

Constant I inches 

0.831 

1.355 

5.200 

(a) Objective function inches 

0.500 

1.000 

11.013 

(b) Weight, lbf 

21.12 

30.24 

212.87 

f r Hz 

1.86 

1.98 

1.46 (e) 

f 2 , Ha 

11.6 (e) 

12.4 

9.17 (e) 

Optimal Solution (c) 




No. of Iterations 

15 

15 

18 

Objective function inches 

0.540 

0.541 

0.547 

Weight lbf 

21.85 

21.87 

21.98 

f jt Hz 

2.00 

1.94 

1.92 

(d) f 2 , Hz 

12.20 

12.20 

12.26 

_*■ , , 4 . inches. 

1^ inches (t^ ) 

1.250(0.0896) 

1.262(0.0908) 

1.327(0.0972) 


(Elements 2 through 10 were essentially all at the lower bound 
value of the side constraint on thickness: t - 0.05 inches, 

I - 0.831 inches^) 


Notes: 

(a) Objective function inches' - £ t. 

i-1 1 

10 

(b) Weight lbf - 12 + 18.24 Z t. 

i-1 * 

(c) All gradients were computed by analytical techniques. 

(d) All solutions have active frequency contraint: f 2 is at its 

lower bound value. 

(e) Initial frequency violates the constraint. 

(f) In this case, a high value of the parameter DABFUN led to premature 
designation of convergence. The solution is nearly optimal, but the 
objective function was still changing in the 3rd significant figure. 

Table 3: Comparison of Solutions Obtained for Various Initial Designs 





*.1 









-48- 

Inltial Values 
Constant t , inches 

4 

Constant I, Inches 
Objective function, inches 
Weight, lbf 
fj_, Hr 
f 2 . Hr 


(a) 


Oblactiva Function 


Using Default Parameters 


0.10 

1.355 

1.000 

30.24 

1.98 

12.4 

(b) 


QRIGINai 

QF P00 « iuAU-* 


Using Tighter Convergence Criteria 


0.54279 



Iteration 10 c 

0.54279 



0.54212 



Iteration 11 

0.54279 

(No 

change) 

0.54076 



Iteration 12 

0.54083 



0.54076 

(No 

change) 

Iteration 13 

0.54033 

(No 

change) 

0.54076 

(No 

change) 

Iteration 14 

0.54083 

(No 

change) 

0.54076 

(No 

change) 

Iteration 15 

0.54072 






Iteration 16 

0.54072 

(No 

change) 




Iteration 17 

0.54072 

(No 

change) 




Optimal Solution 




15 



No. of iterations 

17 



0.54076 



Objective function inches 0.54072 


21.37 



Weight, lbf 

31.87 



1.94 



Hr 

1.96 




12.20 

1.2619(0.09076) 

Notes: 

(a) Default values: 


f 2 . Hr 

T , . 4, inches. 

I^inches (t^ ) 


12.20 

1.2615(0.09073) 


ITEM - 3 
DELFUN - 0.0001 
DAB FUN - 0.001 

(b) Tighter Criteria: ITRM - 5 

DELFUN - 0.00005 
DAB FUN - 0.0005 

(c) The first ten iterations were identical. 


Table 4: Convergence History tot Two Different Convergence Criteria 
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(a) 

9 

No. of 
Iterations 

Final Objective 
Function 

4 

1^ inches 
(■clinches) 

t v Hr (b) 

1 

15 

0.540 

1.262(0.0908) 

1.94 

10 

19 

0.546 

1.315(0.0960) 

2.01 

100 

13 

0.552 

1.308(0.0953) 

2.03 

700 

Kotos: 

26 

0.549 

1.190(0.0838) 

2.06 


(a) For all runs, Default Parameters were used: 

PHI - 5.0 
ITEM - 3 
DELFUN - 0.0001 
DABFUN - 0.001 

(b) The second frequency, f^, was at the lower bound constraint 
(f 2 - 12.2 Hr). 

(c) The initial design was based on constant thickness, constant moment 

4 

of inertia of 0.10 inches and 1.355 inches respectively. 


Table 5: The Effect of Parameter THETA on the Ootimal Solution 
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Initial Value3 


Constant t, inches : 0.25 

Constant I, inches^ : 2.675 

(a) .Objective function jinchcs : 2.50 

(b) Weight ,lbf : 57.60 

f r H = : 2.02 

f 2 , Hz : 12.6 


ORIGIN/"' - 
OF POOR 


p 

QUAi.1T/ 


Optical Solutions (c) 



Case 4 1 

Case 5 f 

Case 6 


f lL " 

1.9 Hz 

f lL“ 

1.5 Hz 

f lL- 

2.1 Hz 


f iu" 

2.1 

f m" 

1.9 

f iu " 

2.5 


f 2 L ” 

12.3 

f 2 L " 

12.8 

f 2 L ” 

11.6 


f 2U " 

12.5 

f 2 U " 

13.2 

f 2U " 

12.0 

No. of iterations 

32 


29 


28 


Objective function, inches 

0.552 

0.662 ] 

0.593 

Weight, lb f 

22.07 

24.08 

22.81 

f r Hz 

2.09 (e) 

1.87 

2.11 (d) 

f 2 » Hz 

12.32 (d) 

12.80 (d) 

11.96 

t^, inches 

0.1022 

0.1122 

0.0753 

t 2 , inches 

0.0500 (f) 

0.0500 (f) 

0.0703 

t^, inches 

/ 


J 


0.0886 

t^ , inches 



\ 

f 

0.0586 

t^ , inches 



0.0500 (f) 

0.0500 (f) 

tg , inches 



0.0623 

J 

. 

t_ inches 
7 1 



0.0791 



tg^ inches 



0.1059 



tg ( inches 

> 

f 

0.0515 


f 

ti 0 , inches 

0.0500 (f) 

0.0514 

j 0.0500 (f) 

Notes: 


10 




10 

(a) Objective function inches - E t, (b) 

Weight 

,lbf - 12 

+ 18.24 

v t 
“ c i 



i -1 




1-1 1 

(c) ITMAX - 40, ITRM - 5 

, DELFUN - 0.0001, 

DAB FUN 

- 0.0025 

for case 5 and, 6 . 


For Case 4, ITEM - 3. 

(d) Active lower bound contraint. (e) Active upper bound constraint, 

(f) Active side constraint. 


» o 

j i 

. 

l i 


i 

1 

i 




Table 6: 


Placement of Frequencies 



ORIGINAL p.'-c:: !3 
0 ?. POOR QUALITY 


Initial Values 

Constant t, inches 

4 

Constant I, inches 

(a) Objective function, inches 

(b) Weight ,lbf 
t v Hz 
f 2 , Hz 
f 3 , Hz 
ITMAX 
ITEM 
DELFUN 
DAB FUN 

f lL’ f lU’ H2 
f 2L ,f 2U* HZ 
f 3L ,f 3U’ H2 

(c) Optimal Solution 
No. of iterations 
Objective function inches 
Weight lbf 

f 1 . Hz 
f 2 , Hz 

f 3 , Hz <d 

t^ inches 

t 2 inches (f 

t^ inches 
t, inches 
tj inches 

tg inches (f 

t 7 inches 
tg inches 

tg inches (f 


Case 7 

Case 8A 

Case 8B 

0.10 

0.10 

0.05 

1.355 

1.355 

0.831 

1.000 

1.000 

0.500 


30.24 

1.98 

12.4 

34.8 

40 

3 

0.0001 
,©o l 

1 . 8 , 2. 2 
12.2,12.6 
34.6,33.0 

20 

0.579 

22.56 

2.02 

12.51 

34.60 

0.1178 

0.0500 


30.24 

1.98 

12.4 

34.8 


0.0500 

0.0501 

0.0608 

0.0500 

0.0500 


57 

0.623 

23.37 

2.11 

12.00 (e) 
33.62 
0.0906 
0.0560 
0.0771 
0.0995 
0.0500 (f) 


21.12 

1.86 

11.6 

32.7 


80. 

5. 

0.0001 
0.0005 
2 . 1 , 2. 5 
11.6,12.0 
33.6,34.0 


44 

0.611 
23.13 
2.10 (d) 
12.00 (e) 
33.60 (d) 
0.0884 
0.0526 
0.0744 
0.0951 
0.0500 (f) 


0.0500 (f) 


0.0500 (f) 


( c o iT-t I rw 
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Notes : 2.0 

(a) Objective function, inches ■ It. 

i»l 

10 

(b) Weight, lbf - 12 + 18.24 Z t 

i-1 1 

(c) All gradients were computed by analytical techniques 

(d) Activa lower bound constraint 

(e) Active upper bound constraint 

(f) Active side constraint 

Table 7A: Optimization with Three Frequency Constraints 









.-.■*** . 


ORIGINAL FA 

G2 t3 

■ • . ■ • ' .V-< 



OF POOR QU 

ality 
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Initial Values 

Case 9A 

Case 9B 

Case 9C 


Constant t. inches 

0.10 

0.10 

0.25 

■ 

* 4 

Constant I, inches 

1.355 

1.355 

2.675 

J 

(a) Objective function, inches 

1.000 

1.000 

2.500 


(b) Weight, lbf 

30.24 

30.24 

57.60 

1 

fj_. Hz 

1.98 

1.98 

2.02 

\ • t 

| 

f 2 , Hz 

12.4 

12.4 

12.6 

1 

1 

f 3 , Hz 

34.8 

34.8 

35.5 

i 

ITMAX 

40 

40 

80 

f ** 

i 

ITEM 

3 

5 

5 

} 

\ 

DEL FUN 

0.0001 

0.0001 

0.0001 

i 

DAB FUN 

0.001 

0.0005 

0.0005 

i 


1L * f 1U ’ 

Hz 

< 

1.5, 1.9 

^ 

2L ,f 2U’ 

Hz 

-< 

12.8,13.2 

* 

3L’ f 3U’ 

Hz 

-4 

35.6,36.0 

* 


(c) Optinal Solution 
No. of iterations 
Objective function inches 
Weight, lbf 


f Hz 
f 2 , Hz 
f 3 , Hz 
t^ ( inches 
t 2 , inches 
t^, inches 
t, .inches 
tj, inches 
tg, inches 
t^ , inches 
tg .inches 
tg .inches 
tlO | inches 


15 

0.697 
24.72 
1.87 (e) 
12.80 (d) 
35.70 
0.0927 
0.C500 (f) 

i 

0.0500 (f) 
0.0637 
0.0827 
0.0917 
0.1064 
0.0602 
0.0500 (f) 


25 

0.684 
24.47 
1.86 
12.83 
35.79 
0.1025 
0.0500 (f) 


'' 


0.0500 

(f) 

0.0705 


0.0943 


0.1166 


0.0500 

(f) 

0.0500 

(f) 


47 

0.684 
24.47 
1.85 
12.82 
35.73 
0.1074 
0.0500 (f) 

t 

0.0500 (f) 
0.0621 
0.0808 
0.1335 
0.0500 (f) 
0,0500 (f) 


( ccm^ i 



t- ■ -v, v^f.'' t 'i-i s 1 ,< .w *■• js' . •■ "-‘-v ' : ■■*'■ .'.. , - : ' -*.vx ^ •*■ :v;-. ■■ >*-»i : >' 1 -Vt."* : ■; " i: •' :% .V 


No tea: 
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10 


(a) Objective function. inches « E t. 

i-1 1 

10 

(b) Weight, lbf - 12 + 18.24 E t. 

i-1 


(c) All gradients ware computed by analtyical techniques 

(d) Active lower bound constraint 

(e) Active upper bound constraint 

(f) Active side constraint 

Table 7B: Optinizatlon with Three Frequency Constraints 
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Initial Values 

Caen 10 


Cana 11A 


Case 11B 1 

Constant t, inches 

< 


- 0.10 


y 


A 

Constant I •inches 

< 


— 1.355 


> 

] 

Density, y^ lb*V inches’* 

0.05 

t 

0.10 


0.10 

j 

Constant W^ ( lbf 

1.512 


1.000 


0.000 

i 

\ • j 

(a)Objectivo function jlbf 

24.24 


28.24 


18.24 

i 

I 

(b)WQijht, lbf 

30.24 


40.24 


30.24 

i 

4 

! 

fj. Hz 

1.98 


1.66 


1.98 

j - ; 

f 2 , Hz 

12.4 


10.7 


12.4 

1 I 

f 3 , Hz • 

34.8 


29.9 


34.8 

| i 

f lL ,f H!' ES 

* 

— 

L.8,2.2 


> 

1 

f 2L ,f 2U’ 33 

< 

— 12.2,12.6 


> 


f 3L ,f 3L” H " 

< 

— 34.6,35.0 


> 


(c)Optimal SoLution 







No. of iterations 

53 


38 


19 


Objective function lbf 

9.252 


10.557 


10.557 


Weight lbf 

15.25 


22.56 


22.56 


f r Hz 

2.12 


1.96 


2.03 


f 2 , Hz 

12.60(a) 


12.51 


12.51 


f,, Hz 

35.00(e) 


34.60(d) 


34.60(d) 

t^inches and W^ lbf 0.0500(f) & 0(f) 

0.1150 & 0(f) 

0.1147 & 0(f) 

tjinches and W 2 lbf i 

^ 0(f) 

0.0500(f) ‘ 


0.0500(f) 


t^inches and W g lbf 

0.056 




i 



t^ Inches and lbf 

2.214 







tj inches and Wj lbf 

0.001 




' 

f 


flinches and W g lbf 

0(f) 

0.0500(f) 


0.0500(f) 


t^inches and lbf 

1.194 

0.0535 


0.0525 


tginchea and W g lbf 

0(f) 

0.0618 


0.0616 


Clinches and W lbf ' 

f 0(f) 

0.0500(f) ' 

' 

0.0500(f) n 

r 

.^inches and W^lbf 0.0500(f) 1.227 

O. 0500(f) «>Q(f ) 

0.0500(f) & 0(f) 


( (T on V I'wc J 
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Notes: 

(a) Objective functionjlbf 


10 

102.4 y t t. + 
i-1 1 


10 
Z W 
i-1 


i 


(b) Weight, lbf - 120y + objective function 

(c) All gradients were computed by analytical techniques 

For all cases, UMAX - 80, ITEM - 5, DELFUN - 0.0001, DAB FUN - 0.001 

(d) Active lower bound constraint 

(e) Active upper bound constraint 

(f) Active side constraint 


Table 8A; Optimization Including Lumped Weights 


Initial Values 


Casa 12 


Constant t,inchaa 

4 

Constant I. inches 

3 

Donaity | Y jlbw/inchao 
Constant W i( lbf 

(a) Objective function, lbf 

(b) Haight, lbf 


0.10 
1.355 
0.10 
1.000 
28. 24 
AO. 24 


f l’ 

He 


1.66 

<r 

Hz 


10.7 

V 

Hz 


29.9 


f lU* 

He 

2. 6, 3.0 

f 2L’ 

f 2U’ 

Hz 

10.5,11.0 

f 3L* 

f 3U' 

Hz 

31.6,32.0 


(c) Optimal Solution 
No. of iterations 
Objective function, lbf 
Haight, lbf 


f l* 

Hz 




i v 

Hz 




f 3* 

Hz 




l l 

Inches 

and 

W l. 

lbf 

c 2 

lnchas 

and 


lbf 

C 3 

inches 

and 

«3. 

lbf 

C 4 

inches 

and 

W 4> 

lbf 

C 5 

lnchas 

and 

W 5' 

lbf 

C 6 

inches 

and 

W 6, 

lbf 

C 7 

inches 

and 

W 7. 

lbf 

C 3 

inches 

and 

V 

lbf 

C Q 

inches 

and 

V 

lbf 

c io 

Inches 

and 

W 10 

lbf 

1 


46 

23.697 
35.70 
2.60 (d) 
10.99 


31.99 


0.1360 & 

0.527 

0. 172o 

0.435 

0.1976 

0.426 

0.2053 

0.339 

0.0769 

1.395 

0.0500 ( f , 

1.135 

\ 

0 


<n 




0.0500 u) 


V 

0 if) 



Notes: 


10 • 10 

(a) Objoctiva function. lbf ■ 182. 4y t t. + E H 

i-1 1 i«*l 1 

(b) Weighty lbf • 120y +• Objective function 

(c) All gradients vera computed by analytical techniques 
ITMAX - 80. ITKM - 5, DELFUH » 0.0001, DAD FUN - 0.001 

(d) Active lover bound constraint 
(o) Active upper bound constraint 
(f) Active side constraint 

Table 8B: Optlniaation Including Lunped Weights 


ORIGINAL PASS 13 
OF POOR QUALITY 


Initial Valua 9 



Case 10 

Case 13 


4 

Constant t inches, I inches , 

3 

Y,l.br/ inches 

0.10, 1.355, 

, 0.05 


Constant W^ lbf. 



1.512 



(a) Objective function, lbf 



24.24 



(b) Weight, lbf 



30.24 



f Hz 



1.98 



f 2 , Hz 



12.4 



f 3 . Hz 



34.8 



f lL’ f 10’ Hz 



1 . 8 , 2. 2 



f 2L* f 2U’ Hz 



12 . 2 , 12.6 



f 3L’ f 3U’ Hz 2 



34.6,35.0 



I^j^lbf inches sec 



O. 

.1100 


(c) Optimal Solution 






No. of iterations 



53. 

80 


Objective function lbf 



9.252 

11.596 


Weight lbf 



15.25 

17.60 


f r Hz 



2.12 

1.80 (d) 


f 2 , Hz 



12.60 (e) 

12.59 


f 3 , Hz 



35.00 (e) 

35.00 (e) 


2 

I lbf inches sec 
m 



807.7 

1144 (g) 


t^ inches and W^ lbf 

0.0500 

(f) * 0 (f) 

0.0500 (f) & 

0 (f) 

t 2 inches and W 2 lbf 

t 


0 (f) 

0.0869 

0 (f) 

tg inches and Wg lbf 



0.056 

0.0757 

0.131 

t^ inches and lbf 



2.214 

0.0500(f) 

0.516 

tj inches and W^ lbf 



0.001 


0.235 

tg inches and lbf 



0 (f) 


0.142 

ty inches and W^ lbf 



1.194 


0.719 

tg inches and Wg lbf 



0 (f) 


1.175 

tg inches and Wg lbf 


/ 

0 (f) 

0.0500(f) 

1.174 

C 1Q inches and W lbf 

0.0500(f) 1.227 

0.1094 

1.833 
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Notes: 

(a) Objective function ( lbf 


10 10 

182. 4y Et.+ S W 
i-1 1 i-1 ' 


(b) Weight lbf - 120y + objective function 

(c) All gradients were coaputed by analytical techniques 

For all cases, ITHAX - 80, ITRM - 5, DELFUN - 0.0001, DABFUN - 0.001 

(d) Active lower bound constraint 

(e) Active upper bound constraint 

(f) Active side constraint 

(g) Very close to being active auto-rotational constraint 


Table 9: Optimization Including Auto-Rotational Constraint 
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0R1G1MAL PAG2 tU 
OF POOR QUALITY 


3 

Y - 0.05 lbs\/ Inches . t ■ 0.10 inches 
(a) W i - 1.5 lbf 

(a) W ± - 100 lbf 


fl ■ 0 RPM 
f x - 2.0 Hz 

t 1 - 0.3 Hz 


ft - 300 RPM 
■ 5.7 Hz 

f x - 5.1 Hz 


3 

Y * 0.05 lbn/lnche3 , t - 1.25 lnche3 
(a) - 1.5 lbf 

(a) W ± - 100 lbf 


ft ■ 0 RPM 
f L - 1.4 Hz 

f x - 0.7 Hz 


fl “ 300 RPM 
f^ “ 5.5 Hz 

f 1 - 5.2 Hz 


Note: (a) Lumped weights were assumed to be uniformly distributed 

for purposes of computation. 


Table 10: Influence of Rotational Speed, Mass, and Stiffness on the 

Fundamental Frequency of a Uniform Cantilever 
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ORIGMjzl I-.'-' r-. 

of poor quality 


Initial Values 


Case 14 


Casa 15 


Case 16 


Constant, t, Inches, I Inches , y lbf/cu.ln 
Conotant W^lbf 

(a) Objoctiva function, lbf. 

(b) Weight, lbf 


•1L’ 


f iu’ Ha 


f 2L’ f 2U’ Q= 

f 3L* f 3U* Hs 2 

I . .lbf inches see 
min' 

n,RPM 


(c) Optinal Solution 
No. of iterations 
Objoctiva function, lbf 
Weight, lbf 
f r Hi 

f 2 » H2 

f v Hr 

J 2 
I . lbf inches sac 
m' 

inches and W^ ( lbf 
t.,, inches and W,,lbf 
t^, inches and W^lbf 
t, .inches and W, .lbf 
tj, inches and W^lbf 
t^, inches and W^, lbf 
t^, inches and W^,lbf 
t„ inches and W_ lbf 

O' 0 » 

tg finches and W^lbf 

t,„ inches and W,- lbf 
lOi lOi 


68 

9.252 
15.25 
2.12 
12.60(e) 
35.00(e) 
796.7 

0.0500(f) & 0 (f) 

0 (f) 

0.022 
2.377 
0 (f) 

0 (f) 

1.145 
0 (f) 

^ 0 (f) 

0.0500(f) 1.149 


A 


0.10,1.355,0.05 
1.512 


— 24.24 - 

- 30.24 - 

- 1 . 8 , 2. 2 • 

12.2,12.6 

34.6,35.0 

500 — 

30 


76 

9.303 

15.30 

2.16 

12.60(e) 

35.00(e) 

811.7 

0.0500(f) 6 0 
* 0 
0 


V 

0.0500(f) 


(f) 

(f) 

(f) 


2.345 
0 (f) 

0 (f) 

1.121 

0 (f) 

0 (f) 

1.277 


100 


80 

16.718 

22.72 

2.20(e) 

12.59 

34.97 

1957 

0.0500(f)&0(f) 


V- 

0(f) 

i.oo: 

0 (f) 
0. 40t 
y 1.44 f 

0.0500(f) 0(f) 
0.3685 6.39? 


(do* iKrtvetY \ 
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Notea : 

(a) Objective function lbf 


10 10 
132. 4V* E t. + E W 
i-1 1 i°l 


(b) Weight lbf “ 120y + objective function 

(c) All gradients vare coaputed by analytical techniques 
Eigenvalues were computed with double precision routines 

For all cases, ITMAX - 80, ITEM - 5, DELFUN - 0.0001, DAB FUN - 0.001 

(d) Activa lower bound constraint 

(e) Active upper bound constraint 

(f) Active side constraint 

(g) Active auto-rotational constraint 


Table 11A: Optimization of Rotating Bean 
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Initial Values 

4 

Constant t , inches , I inches , ylb«/cu. in. 
Constant W^,lbf 

(a) Objective function j Ifa-f. 

(b) Weight f lbf 

f lL» f lu’ 02 

f 2L’ f 2U’ Hz 
f 3L , f 3(J> Hz 

^min ^ b ^ Inches sec 

n RPM 

(c) Optinal Solution 
Wo. of iterations 
Objective function lbf 
Weight lbf 
f 

f 2 , Hz 
f 3 , Hz 

I lbf inches 

tQ 

t^, inches and 
t 2 , inches and 
t 3 , inches and 
t^ , inches and 
tj, inches and 
t... inches and 
t^, inches and 
tg, inches and 
tg, inches and 
t^Q, inches and 


Case 17 

0.10,1.355,0.05 

1.512 

24.24 

30.24 
5. 2,5. 6 

18. 0,18. 6 
42.0,43.0 
1100 
300 


39. 

10.539 
16.54 
5.59 (e) 
18.59 (e) 
42.99 (e) 
1113. (g) 


0.0500 

(f> & 

0 

(f) 

1 



0 

(f) 




0.002 





0.591 





0.712 





0 

(f) 




0.745 





0.438 


\ 



0 

(f) 

0.0500 

(f> 

3.492 



sec 

W^lbf 

w 2 , lbf 

W 3 ,lbf 
W. .lbf 
U 5 \lbf 

W 6, lbf 

W 7( lbf 

W 8 . ^ 
Wg.lbf 

w iof bf 


[covdt fvwe^N 
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Notes: 

(a) Objective function jlbf 


10 10 
182. 4y Z t. + Z W. 
i-1 1 i-1 1 


(b) Weightj lbf «■ 120y + objective function 

(c) All gradients were computed by analytical techniques 
Eigenvalues were computed with double precision routines 

For all cases, ITHAX - 80, ITEM - 3, DELFUN - 0.0001, DAB FUN - 0.01 

(d) Active lower bound constraint 

(e) Active upper bound constraint 

(f) Active side constraint 

(g) Active auto-rotational constraint 


Table 11B: Optimization of Rotating Beam 
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Frequency 


No. of 
Elements 

First 

Mode 

(Hz) 

Error 

. _(*> . 

Second 

Mode 

(Hz) 

Error 

(X) 

6 

1.9678 

0.6 

12.142 

3.6 

8 

1.9733 

0.3 

12.258 

2.7 

10 

1.9761 

0.2 

12.312 

2.3 

12 

1.9780 

0.1 

12.342 

2.0 

14 

1.97C9 

0.06 

12.360 

1.9 

20 

1.9805 

0.03 

12.385 

1.7 

Analytical 

1.980 


12.6 



Table 12; Frequency vs. Number of Elements 
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Frequency 

(c?e) 


Mode No. 

(1) 

(2^ 

(3) 

(4) 

(5) 

1 

2.01622 

1.95462 

2.01614 

2.01614 

2.01511 

2 

12.6354 

12.6334 

12.6355 



3 

35.3872 

35.3916 


35.3875 


4 

69.3934 



69.3937 


5 

114.892 



114.892 


6 

172.121 

Rest are 

essentially 

172.122 

rest are 

7 

241.495 

Identical 

. to (1) 

241.495 

essentially 

8 

323.537 



323.536 

identical 

9 

418.329 



418.329 


10 

520.042 



520.043 

to (1) 

11 

692.152 



692.155 


12 

836.210 



836.210 


13 

1013.68 



1013.63 


14 

1223.45 



1223.44 


15 

1470.15 



1470.15 


16 

1758.30 



1758.31 


17 

2086.84 



2086.85 


18 

2436.29 



2436.30 


19 

2743.53 



2743.33 


20 

3433.47 



3433.43 



Table 13: Frequencies Calculated bv Various Methods 




\ 
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Eleraent Thicknoaaas 
(inches) 


Element 

No. 

Single-Preclaion 

Calculation 

Double-Precision 

Calculation 

1 

0,0753 

0.0724 

2 

0.0703 

0.0881 

3 

0.0836 

0.0821 

4 

0.0856 

0.0564 

5 

0.0500 

0.0500 

6 

0.0500 

0.0500 

7 

0.0500 

0.0500 

8 

0.0500 

0.0500 

9 

0.0500 

0.0500 

10 

0.0500 

0.0500 

Objectivo 

Function 

(LBF) 

0.593 

0.599 


Table 14; Difference In Optimal Designs Caused bv Errors In 


Eigenvalue Calculation 




i 

i 

i 


-69- 


Elcnant 

Ho. 

(1) 

3 V 3t i 

(cps/ln.) 

(2) 

3f 1 /3Q l 

(cpa/lbf) 

(3) 

(cpa/lbf) 

1 

4.791 

0.525 

0.000 

2 

3.492 

0.383 

-0.001 

3 

2.388 

0.262 

-0.003 

4 

1.472 

0.161 

-0.009 

5 

0.729 

0.080 

-0.022 

6 

0.105 

0.012 

-0.043 

7 

-0.443 

-0.049 

-0.073 

8 

-0.957 

-0.105 

-0.114 

9 

-1.489 

-0.163 

-0.166 

10 

-2.074 

-0.227 

-0.229 

Table 15: 

Sensitivitv of Frequency 

to Changes in Thickness 

and Weight 



II, FIGURES 
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originai pr ey, a 

OF POOR QUALIP,' 



Figura 2: Element undor Constant Tension. 



TOTAL NUMBER OF ELEMENTS IN MESH 


Figure 3: Obj vs. No. of Elements. 



WEIGHT (LBF) 
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ORIGiNAL IS 
Op POOR QUALITY 



n, TOTAL NUMBER OF ELEMENTS IN MESH 


Figure 4: Lumped Weight Versus No. of Elements. 



n, TOTAL NUMBER OF ELEMENTS IN MESH 


Figure 5: Moment of Inertia at Free End Versus No. of Elements in Mesh. 



WEIGHT (LBF) 
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0RIG5NAL PAGE 13 
• OR POOR QUALITY 



tTLjS I.8SS0 


(130 


h 


X 




rums 


.m 3 3 5.7746 


OSJ a * 21.6110 
(Noto:08J 3 >OBJ s ) 


-I a = 3.3512 

10.4S03 



Figure 7: Optimal Designs for Various Values of n. 
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ORIGJNAL PAGZ 13 
OF POOR QUALITY 



mg’ 0.92519 [V 


n=l4 m a = 0.78326 


,rr.| 4 = 5.1907 

0BJ i4 = 20.8154 



X 14 = 4.4000 
|j m l4 *m 14 = 9.3064 


Figure 


7: (Continued) 


CRIGlFiAL FAGS E3 
OF FOGS QUALITY 




m 2o = 7.98IS 


08J 2 o= 20.3763 


m 20 + rn 20 = 8-7980 
I 20 = 1.488! 


Figure 7: (Continued) 
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ORIGHVAL PAGc .'S 
OF POOR QUALITY 



Figure 8: 


Alternative Optinum - Found by Imposing Constraint 0.83073 < I 10 < 2-0 
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